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Abstract By extending the minimum spanning tree (MST) clusterization algorithm for the binding energy cut, the 
isospin asymmetry dependence of directed flow for isospin sensitive isobar pairs (neutrons-protons, ?H-^He) is studied 
from low towards high incident energies. The modified clusterization method (MSTB) has the advantage to identify 
the fragments at quite early time. It enhances (reduces) the production of free nucleons (fragments) over MST method. 
The directed flow of isobaric pair *H-*He is more sensitive towards isospin asymmetry caused by MSTB than isobaric 
pair n-p. This sensitivity becomes quite strong towards the high incident energy and neutron-rich reaction system. In 


conclusion, the inclusion of binding energy in clusterization method for the flow studies has been uniquely important 


for understanding the isospin physics, especially for high density behavior of symmetry energy. 
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1 Introduction 


With the availability of radioactive beam facilities, 
isospin physics has attracted the whole of the nuclear 
community in last decade. In intermediate energy 
heavy-ion collisions, the dynamical models such as 
Quantum Molecular Dynamical (QMD)'- Pl, 
Boltzmann-Uehling-Uhlenbeck (BUU)! and 
statistical models^??" were used to get the rich 
information about isospin physics. The dynamical 
models are of unique importance due to their ability 
to follow the time evolution from initialization 
through compression and expansion to equilibrium 
state. It is also well known that dynamical models do 
not simulate the fragments. In order to simulate the 
fragments, one needs the secondary algorithm for 
coupling with dynamical models. In the efforts to 
reproduce the experimental data, many secondary 
algorithms have been developed to study the nuclei 


[11-13,15,28,29 


near the drip line l but very few for nuclei 


< qno [5,6,8,10,15 
away from drip ine 


The most commonly and widely used 
algorithm depends on the spatial and momentum 
coordinates of the nucleons, is known as minimum 
spanning tree (MST) algorithm ^ ^l. According to this 
method, two nucleons undergo the cluster formation if 
the relative distance (|R;-Rj]) and relative momentum 
(|P—P;|) between the nucleons was less than 3.5—4 fm 
and 250—268 MeV/c, respectively. These parameters 
can be obtained by fitting the experimental data for 
some of the global observables such as multiplicity of 


[2,55] with theoretical 


intermediate mass fragments 
results. Recently, MST method was further extended to 
isospin-dependent MST, in which the cut on the 
momentum space was kept same, but the cut on the 
spatial coordinates was constraint on the basis of type 
of particles. The distance between the different kind of 
particles was taken as follow: |R? — R?|— 3 fm, 
|R? — R? |H R’ - R? |= 6 fm", 

Since the MST method only depends on the 
constraints from position and momentum, it seems to 
be worried about the stability of fragments due to the 
formation of artificial weakly bound fragments. To 
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avoid from this problem, more complicated methods 
like Stimulated Annealing Clusterization Algorithm 
(SACA)? Early Cluster Recognization Algorithm 
(ECRA)P?l were also developed. These methods were 
found quite successful, but, very complicated. 
Moreover, due to the choice of parameters such as 
cooling parameters, iterations procedure, and choice of 
minimal can follow the totally different configuration 
of fragments. Due to the limitation of sharp minimal, 
in the mildly excited or asymmetric systems, the scope 
of these two methods was found to restrict at some 
points. The possible best method for avoiding from the 
artificial formation of fragments is found a way to 
constrain the fragments by using an average binding 
energy cut of 4 MeV/nucleonP?! One can further 
improve the method by using the realistic binding 
energy rather than 4 MeV/nucleon". This method 
was found to be as simple as MST and found to 
reproduce the experimental data just like the 
complicated methods SACA and ECRA etc. The 
extensive studies for the multiplicity of different kind 
of fragments with different clusterization methods 


dU 2-152850 


have been performe | while flow parameters 


are poorly known with this type of 
methods U^ 1:1231.32] 

In last two decades, the collective flow of 
nucleons with non-central collisions was found as a 
useful tool for extracting the nuclear equation of state 

(NEOS) of symmetric nuclear matter? ?! With the 
passage of time, efforts have been done to extract the 
information of isospin physics in term of symmetry 
energy using the directed as well as elliptic flgy Ple 
The studies showed that the experimental results of 
directed flow can be reproduced by (1) changing the 
mean field or nucleon-nucleon (NN) cross sections in 
the transport model (2) using the nucleon phase space 


[31334] Most of the studies were found 
33-35 


for calculations 


to obey the hard equation of state"? which opposed 


the findings with other observables in the 


12,19,28,31,32,36-38] 


literature! . In recent study, the soft 


[39] but, the comparison 


equation of state was predicted 
was made between the directed flow extracted from 
phase space of nucleons and the directed flow of 
Within the BUU 


, the importance of directed flow of 


fragments having Z < 2%, 


calculations’! 


isospin sensitive fragments n, p, ?H and *He over the 
directed flow of nucleon phase space has already been 
put forwarded. The elliptic flow extracted with the 
simple MST method was also found to be quite 
sensitive towards the determination of symmetry 
energy at supra-saturation densities"! The results 
indicated that the directed as well as elliptic flow can 
be a useful tool to understand the isospin physics or 
symmetry energy. 

The main problem with the study was that 
isospin physics with directed and elliptic flow is either 
studied from the nucleon phase space or by using the 
simple MST method. It is worth mentioning here that 
one cannot guarantee about the stability of fragments 
formed with simple MST method. Due to the worry 
about the stability of fragments, it is the prime need of 
the present time to study the importance of Isospin 
physics from directed and elliptic flow with the stable 
fragments. 

In the present study, first the stable fragments 
are formed by applying the binding energy cut on the 
pre-clusters formed with the simple MST method. 
Afterward, the sensitivity of the isospin sensitive 
fragments n, p, ^H and ?He directed flow is studied 
from low towards the high incident energy with the 
variation in isospin asymmetry of the reaction systems. 

The modified 
explained in Section-2, the results and discussions are 


clusterization method is 


presented in Section-3, followed by the conclusion in 
Section-4. 


2 Methodology 


In the present work, IOMD model is used, which is 


i s "m s " 5.6.10 
discussed in detail in our recent publications! m 


originally developed by Hartnack and Co-workers!! Pl. 
The model is modified by the authors for the density 


dependence of symmetry energy, having form: 


ek (2 j^ m ( p yt 
2 p 2 p 


Es (p) = 


with the parameters of C,,= 25 MeV and C,,— 35.2 
MeV. When we set 7;=0.5 and 1.5, respectively, it 
corresponds to the soft and stiff symmetry energy" !. 


Minimum  Spanning Tree method with 
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Binding Energy Check (MSTB) is used. This method 
is a modified version of the normal MST and old 
MSTB method. The difference between old MSTB and 
this MSTB is the inclusion of energy from momentum 
dependent interactions as well as symmetry energy 
along with Skyrme interactions. The procedure is as 
follow: the phase space obtained from IQMD is 
analyzed with simple MST method and pre-clusters 
are sort out. Since we are not aware about the stability 
of pre-clusters formed at this stage, the pre-clusters 
formed from the simple MST are now subjected to the 
binding energy condition as follow: 


N 
P, Y +m, =m, ZA <—Epging: 

Bea 
Here, we take Egina = 4.0 MeV/nucleon if Nr23 and 
Egpiga— 0 otherwise. In this equation, Nr is the number 
of nucleons in a fragment, Py, is the average 
momentum of the nucleons bound in the fragment. The 
requirement of a minimum binding energy excludes 
loosely bound fragments which will decay later. The 
realistic value of Egina changes slightly the fragment 
multiplicity at intermediate times, but has no influence 
on the qualitative behavior and on the asymptotic 
results. However, if by using the realistic binding 
energy, one searches for the most bound configuration, 
the results will get affected. At the present time, we 
had just focused on the bound configurations and 
hence the average cut 
—4 MeV/nucleon is justified. 


of binding energy 


3 Results and discussions 


Several thousands of event for the reactions of 
?Snt!?Sn and '“Sn+!4Sn between the incident 
energy 50 and 600 MeV/nucleon for the nearly central 
collisions using the IQMD model coupled with MST 
and MSTB algorithms are simulated. The soft 
momentum dependent equation of state with soft 
symmetry energy and isospin-energy dependent cross 
sections is employed. 

There are two methods in the literature used to 
calculate the directed flow. In the first case, the 
directed flow is extracted from the mid-rapidity slope 
of the P,/A versus rapidity distribution Y.m/Ybeam plots. 
The rapidity distribution is calculated as follow: 


da EG) ep.) 
2 E()-p,G) 
where E(i) and pi) are the total energy and 


Y(i) 


longitudinal momentum of ith particle, respectively. 
The second method is to study the directed 
transverse in-plane flow P? is as follow: 


(p*)- 5 sien YOO, 


where Y(i) is the rapidity distribution as discussed 
above and P,(i) is the transverse momentum of the ith 
particle in x-direction. This P^ is defined over 
entire rapidity region and therefore expected to present 
an easier way of measuring the in-plane flow rather 


than complicated P,/A plots. 


- Yi 20.5 
^ 
ı bz2fm 
A 50 MeV/nucleon 
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Fig.1 Time evolution of different kind of fragments: free 
nucleons (upper), LCPs (middle) and IMFs (bottom) for central 
collisions at 50 MeV/nucleon with MST and MSTB algorithms. 
The left (right) panels are for neutron-poor (neutron-rich) 
'Psnt!?Sn C'4Snt!4Sn) | reaction systems. All the results are 
with the soft symmetry energy. 

In order to check the stability of fragments 
and effect of isospin physics, in Fig.l, the time 
evolution of free nucleons, light charged particles 
(IMFs) 
for neutron-poor and neutron-rich reaction systems at 
50 MeV/nucleon with MST and MSTB algorithms is 
displayed. In the high density phase, the MSTB 


algorithm does not find any fragments with reasonable 


(LCPs) and intermediate mass fragments 


binding energy and hence most of the particles are free, 
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while, there were a lot of artificial fragments with 
MST. After the high density phase is over, it starts 
recognizing the fragments, which are real, bound and 
stable. In all the cases, MSTB helps to identify the 
fragments quite early. From the figure, it is clear that 
MSTB enhances the production of free particles and 
reduces the production of LCPs and IMFs. Moreover, 
with increase in the size of the fragment, MST takes 
less time to match with the MSTB results and also the 
difference between MST and MSTB results goes on 
decreasing throughout the time evolution. 

When the 


neutron-poor and neutron-rich reaction systems, it is 


results are compared for 
found that enhanced production of fragments takes 
place for the more neutron-rich system. The enhanced 
production is quite effective in free nucleons and LCPs 
compared to IMFs. It is also shown by us recently that 
free nucleons and LCPs are more sensitive towards the 
symmetry energy compared to IMFs!?l. In this regard, 
the sensitivity of free nucleons and LCPs with MSTB 
towards the neutron-rich system can also play an 
important role in the prediction of symmetry energy, 
which is a topic of separate discussion. After assessing 
the sensitivity, for the further study, the isospin 
sensitive free nucleons and LCPs, namely, n, p, 3H and 


3 
He are used. 


P/A (MeV/c) 
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Fig.2 Importance of MSTB in flow calculations by plotting 
the rapidity distribution dependence of P,/A for neutrons, 
protons, 3H particles with MST and MSTB algorithms at 
100 MeV/nucleon for '4Sn+'™4Sn reaction system. 

Using the first method, we have plotted in 
Fig.2, the rapidity distribution of P,/A for neutrons, 
protons and 7H particles at 100 MeV/nucleon for the 


us 124 
central collisions of 


Sn reaction system. The 
mid-rapidity slope of P,/A represents the directed flow. 
The flow is found to be affected only for the heavier 
fragments, neutrons and protons flow is least affected. 
This is due to the isotropic (anisotropic) distribution of 


free nucleons (LCPs)?l. Naturally, yield has to be 


independent of the type of distribution, while 
momentum of particles is supposed to depend strongly 
on the type of distribution. One of the examples is that 
different behavior of flow is obtained in earlier studies 


4241] The sensitivity 


for different distribution regions! 
of flow with MSTB towards the isospin physics at 
higher incident energy is presented in the following 
sections. 
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Fig.3 (Color online) Isospin asymmetry dependence of 
directed flow for isobaric pair neutrons, protons with MST and 
MSTB method. The different panels are at different incident 
energy ranging from 50 to 600 MeV/nucleon. 

In order to study the role of MSTB at higher 
incident energies and towards isospin physics, in Fig.3, 
the isospin asymmetry dependence of more reliable 
flow quantity P^' is plotted for isobaric pairs n, p 
CH He) with MST and MSTB algorithm from 50 to 
600 MeV/nucleon. In Fig.3, with increase in the 
incident energy, isospin asymmetry dependence of 
directed flow for n, p is found to affect more by the 
method of clusterization in comparison with the type 
of particles. In all panels, protons (neutrons) have 
more (less) positive directed flow with MST as well as 
MSTB method. This is due to the contribution of 
Coulomb repulsion in protons production. The directed 
flow of protons as well as neutrons becomes more 
positive with MSTB over MST towards the high 
incident energy. This is true because there was 
increased content of unstable fragments with MST 
towards the high incident energy. The MSTB will 
break these unstable fragments into free particles (as 
per the condition specified earlier) and hence 
dominance of collisions will take place. Moreover, 
higher positive flow of neutrons towards the high 
incident energy can act as a tool for isospin dependent 
cross section. Furthermore, these findings are 
suggesting a shift in the balance energy with different 


kind of fragments. 
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(a) 50 MeV/nucleon 


(b) 100 MeV/nucleo 


— "M 


MST MSTB 
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NIZ 
Fig.4 (Color online) Same as ETES but for the ?H and *He 
particles. 

In Fig.4, just like n and p, the isospin 
asymmetry dependence of directed flow for ^H, ?He is 
found to be affected more by the method of 
clusterization compared to type of particles towards 
the high incident energy. Up to 100 MeV/nucleon, *He 
CH) particles have more (less) positive directed flow, 
but, at and after 200 MeV/nucleon, *H C He) particles 
have more (less) positive flow with MST as well as 
MSTB methods. This is due to the dominance of mean 
field (collisions) at and below (at and above) 100 (200) 
MeV/nucleon. When the mean field (collisions) 
dominates, Coulomb interactions (isospin dependent 
cross sections) make the flow more positive for proton 
rich (neutron-rich) ?He (H) particles. Towards the 
higher incident energies, the flow of ^H as well as *He 
becomes more positive with MST over MSTB method. 
As with the inclusion of binding energy restrictions in 
clusterization algorithm, most of the 3H and ?He were 
found unstable. Especially, this is true for ?He and 
hence directed flow is decreased with MSTB, which 
has resulted in the increase of the transverse directed 
flow for neutrons and protons (as discussed before in 
Fig.3). 

In brief, the directed flow of isospin sensitive 
particles n, p, 3H, ?^He is affected by the different 
method of clusterization. The role of different 
clusterzation method is dominating for the isobaric 
pair CH, *He) compared to (n, p). Moreover, the 
behavior of directed flow for (n, p) pair is smooth with 
incident energy, but, the transition is observed for the 
directed flow of (^H, *He) pair with incident energy. 
The isobaric pair (n, p) and CH, *He) has opposite 
sensitivity with the type of method of clusterization. 


These points suggest that the methods of clusterization 
are going to play an important role in understanding 
the isospin physics and hence need to be handled 
carefully during the determination of symmetry energy 


also. 


4 Conclusion 


The modified clusterization method has a time saving 
advantage in terms of providing the real, bound and 
stable fragments at quite early time. It has the 
capability to obtain the results similar to the 
complicated secondary algorithms SACA and ECRA. 
have the 
drawback over MSTB of time consuming and consists 


Moreover, the complicated methods 
a lot of adjustable parameters. The directed flow of 
isobaric pair ?H-He in comparison to n-p is found to 
be strongly dependent on the isospin asymmetry 
caused by MSTB with the increase of incident energy 
as well as isospin of the reaction systems. This 
indicates it as a robust approach for understanding the 
high density behavior of isospin physics or symmetry 
energy. 
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